BLOWUP OF CLASSICAL SOLUTIONS FOR A CLASS OF 3-D 
QUASILINEAR WAVE EQUATIONS WITH SMALL INITIAL DATA 
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Abstract. This paper is concerned with the small smooth data problem for the 3-D nonlinear wave equation 

3 

d^u — (1 + u + dtu) Au = 0. This equation is prototypical of the more general equation y * g$j (u, ^7u)d-ijU = 0, 

i,3 = 

3 

where xq = t and gij(u, Vit) = cij +diju + y ^ e^jdku + 0(\u\ 2 4- |Vu| 2 ) are smooth functions of their arguments, 

3 

with Cij,dij and being constants, and dij ^ for some (i,j); moreover, y * e^{d^u)diju does not fulfil 

i,j,k=0 

the null condition. For the 3-D nonlinear wave equations d^u — (1 + u) Am = and dfu — (1 + dtu) Au = 0, 
H. Lindblad, S. Alinhac, and F. John proved and disproved, respectively, the global existence of small smooth 
data solutions. For radial initial data, we show that the small smooth data solution of d^u — (1 4- u + dtu) Au = 
blows up in finite time. The explicit expression of the asymptotic lifespan T £ as e — ¥ 0^ is also given. 



§1. Introduction and main results 
We consider the second-order nonlinear wave equation in [0, oo) x W 1 

n 

X! gij(u,S7u)dijU = 0, 

i,j=o (1-1) 
(u(0, x),d t u(0, x)) = (eu (x),eui(x)) , 

where xq = t, x = (xi, ...,x n ), V = (do, 9i, d n ), e > is a sufficiently small constant, uo(x), u\(x) £ 
C oo (M n ), and the g-ij{u, Vw) are smooth functions of there arguments which are of the form 

n 

9ij(u,Vu) = Cij + d ij u + J2e* j d k u + 0{\u\ 2 + |Vu| 2 ) (1.2) 

fe=0 

n 

with cy, dij and ejj- being constants. We assume that the linear part Cijdiju is strictly hyperbolic with 

i,j=o 

respect to time t. From [8-9, 14-16] we have that (1.1) has a global smooth solution when n > 4. 
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If n = 3 and dij = for all < i, j < 3 in (1.2), then (1.1) has a global smooth solution if the null 

3 

condition for the the main part ^ e^dkudijU holds, otherwise the solution of (1.1) blows up in finite 

i,j,k=0 

time. See [4, 8-13, 17, 20-27] and the references therein. 

If n = 3 and dij ^ for some but e\j = for all < i, j, k < 3 in (1.2), then it follows from the 

results in [3] and [18-19] that (1.1) has a global smooth solution. 

3 

The following interesting problem naturally arises: If n = 3, dij ^ for some (i, j), and e^dkudijU 

in (1.2) does not fulfill the null condition, does the smooth solution of (1.1) blow up in finite time or does it 
exist globally? In this paper, we are concerned with this problem, especially (and without loss of generality) 
the protypical equation d^u — (1 + u + d t u)Au = is studied. More specifically, we consider the problem 

| dfu-(l + u + d t u)Au = 0, (t, x) e [0, oo) x R 3 , 
\ (u(0, x),d t u(0, x)) = (suo(x),eui(x)), 

where uo(x), U\(x) € C^°(K 3 ) are radially symmetric and suppuo U suppui C {x : \x\ < M} with M > 
a constant. For notational convenience, we write (uo(r), ui(r)) instead of (uq(x),ui(x)) later on and the 
domains of definition of «o( r ) an d u±(r) are simultaneously extended to [— M, M\. This results from the 
fact that uo(r) and ui(r) are actually smooth functions of r 2 due to uq(x),ui(x) G C^R 3 ). 

Let Fq(s) = - (^suq(s) + J sui(s) ds^j for s£i According to Theorem 6.2.2 and (6.2.12) of [9], we 

know that the function F (s) ^ unless both Uo(s) = and Ui(s) = 0. Moreover, F (s) = for \s\ > M. 

Let r { s) = ^ F ^% s) *» - € M], F>(s) + 0, and > «• Further let 

A = € (— M, M) : F*(s) ^ 0, pll ^ ^ > 0,t(s) > oj and B = { S e (— M, M) : f» = 
0,^' («) > 0}, and denote 

r = min< minr(s), min —— — \. (1.4) 

\seA v h seB F£{s) J v ; 

It can be shown that To is a finite positive number if (uo(r) , U\(r)) ^ holds. 
The main result of this paper is: 

Theorem 1.1. Assume that Uo(x), Ui(x) € C^K 3 ) only depend on r = \/x\ + x\ + x\. If Uq(x) ^ or 
u\(x) ^ 0, i/ien problem (1.3) has a C°° solution u(t, x) for <t < T e , where T e stands for the lifespan of 
the smooth solution u(t,x), which satisfies 

limelnT e = t > 0. (1.5) 

e— >0 

Remark 1.1. It follows from Theorem 1.1 that the smooth solution of (1.3) blows up in finite time provided 
that («o(r), ui(r)) ^ 0. 

Remark 1.2. We have asserted that To > is a finite number as long as (uo(r) , u\(r)) ^ 0. To prove this, 

2 

it suffices to show that T\ = min t(s) > is finite, since T2 = min — — — > obviously holds. We first 

seA seB Fq(s) 

show A ^ 0. Let A 1 = {s e (— M, M) : Ftf(s) > 0}. Then obviously A C A x . Since A x is a bounded, 

open, and nonempty set by (u (r) , Ui(r)) ^ and the continuity of Fq(s), one can write Ai = ^J(a;,fy), 

l 

where the union is finite or countable infinite and {ai 1 , bij and (a; 2 , 6; 2 ) are disjoint for /i 7^ ?2- Moreover, 
F^'(ai) = F£(bi) = for any I. Set B = {s e [-M,M] : F^(s) = 0}. Then B is closed and (oj,6j) \ B 

2 



is an open and nonempty set. Write (ai,bi) \ Bo = ^J(a; m ,6; TO ), where F^aim) = Fg(6; TO ) = and the 

rn 

intervals (a; m ,6; m ) are disjoint for different (l,m). For any (7 ,TOo), there exists at least exists one point 
so e {ai omo ,bi omo ) such that F^'(sq) - F^(s ) > 0. If not, then F^'(s) - F^(s) < for s e (o; omo , bi omo ), 
which means (e~ s Fg(s))' < and then Fq(s) = because of F^aim) = F^bim) = 0. This contradicts 
F£(s) > for s e (ai mni^omj ^~ A\. On the other hand, t(sq) > obviously holds, and thus the set 

A is nonempty. Next we show n > 0. Set z = — jSrrl for s e A. Then r(s) = ^ — — with 

F (s) F (s) z 

s e A, z > -1 and z^O. If z £ (-1, 0), then r(s) > _,„ , , . If z e (0, N] for some fixed > 0, then 

maxiv (s) 

s£A 

Cn 

t{s) > „ , ■ If z > is large, then Fq'(s) > is small and -Fq( s ) < holds. In this case, we conclude 

max _tq ys J 

2 

from t(s) = — ln(l + z) that, if |-Fo( s )l ^ for some fixed constant 5q > 0, then r(s) > C((5o) > 0; 

on the other hand, if |-Fo( s )l i s small, then r(s) > is large. Consequently, putting everything together, 
the constants T\ and To are positive and finite. 

Remark 1.3. The result of this paper can be extended to the more general nonlinear wave equation 

3 

9ij(u,Vu)d ijU = 0, 
i,j=0 ( ' i 

(u(0,x),d t u(0,x)) = (euo(x),eui(x)), 

where xq = t, Uq(x), Ui(x) € Cq°(M™), the gij(u, Vit) are smooth functions of their arguments which are of 
the form 

3 

gij {u, Vu) = + dijU + ^2 e iAu + 0(\u\ 2 + \Vu\ 2 ) 

k=Q 



3 

r,2 



with Cy, <ijj and e^- being constants, and 7^ for some (j, j), and the null condition for e^dkud, 

i,j,k=0 

does not hold. Because a proof of this statement just requires the methods and estimates used in this 
paper and the ones of [5] and [25] combined with blowup system techniques of [1-2], but it is technical and 
tendious otherwise, it is omitted here. 

Remark 1.4. For the 2-D case of problem (1.1), if the coefficients gij(u, Vu) are independent of u, then 
there is a rather complete collection of results on the global existence and the blowup, respectively, of small 
smooth data solutions, see [1-2, 7] and the references therein. On the contrary, if the coefficients gij(u, Vu) 
depend on both u and Vw, there have been no systematic studies so far. Related results will appear in a 
forthcoming paper of ours. 

In order to prove Theorem 1.1, we first derive the lower bound on the lifespan T e for problem (1.3) when 
the initial data are radial. By constructing an approximate solution as in [9] or [6], then considering the 
difference of the exact solution and the approximate solution, applying the Klainerman-Sobolev inequality, 
and establishing some further energy estimates, we obtain a lower bound on the lifespan T e . Here we point 
out that although elnT e > C > has already been shown in [17], for the reader's convenience and to 
obtain the sharp lower bound t , we still give a complete proof. On the other hand, it follows from radial 
symmetry of the initial data (uo(x), ui(x)) that the solution u(t, x) is also radially symmetric for t < T e . 
Based on this, we change (1.3) into a 2 x 2 system of two independent variables (£, r). Then, by using the 
properties of the approximate solution constructed above and some delicate analysis, and by treating the 
solution u accordingly, we obtain the upper bound on T s . Here the derivation is motivated by the methods 
of [10], where the equation d 2 u — c 2 (dtu)Au = was studied, the coefficient c 2 (dtu) of which only depends 

3 



on the gradient of the solution u. Compared with [10] and [6], due to the simultaneous appearance of u 
and dtu in the coefficients of equation (1.3), we have to introduce a few more quantities in order to get a 
"blowup" -type nonlinear second-order ordinary differential equation with suitable initial data that provides 
the upper bound on ehiT e . Based on the results in the two steps above, we finally complete the proof of 
Theorem 1.1. 

In this paper, we will use the following notation: 

Z denotes one of the Klainerman vector fields in the radially symmetric case, i.e., 

d r , d t , S = td t + rd r , H = rd t + td r , 
d stand for d r or d t , and the norm ||/||z,2 means \\f(t, -)IU 2 (K 3 )- 

§2. The lower bound on the lifespan T £ 

In this section, we establish the lower bound on T e for smooth solutions of the Cauchy problem (1.3). As 
in the proof of [9, Theorem 6.5.3], by constructing an approximate solution u a of (1.3) and then estimating 
the difference of u a and the solution u, we obtain the lower bound on T e by a continuity induction argument. 
The new ingredients in this procedure are the construction of the approximate solution and treating the 
solution u itself that occurs in the equation in (1.3) rather than the derivatives of this solution only as in 
[9]. Although this procedure is analogous to the one in [6], for the reader's convenience and also as it is 
applied to obtain the upper bound on T e , we still give a complete proof. 

Let the slow time variable be r = eln(l + t), and assume the solution of (1.3) is approximated by 

-V(q,r), r>0, 
r 

where q = r — t and V(q, r) solves the equation 

2d qT V + Vd 2 q V - 8 q Vd 2 q V = 0, {q, t)£Rx [0, oo), 

V(q 7 0) = F o (q), (2.1) 
, supp V C{q< M}. 

Lemma 2.1. (2.1) has a C°° solution for <t < tq, where To is given by (1.4). 
Proof. Set w(q,r) = d q V(q,r). Then it follows from (2.1) that 

| 2d T w + (V - w)d q w = 0, (q,T)e (-oo, M] x [0,r ), 
The characteristics of (2.2) starting at the point (0, s) is defined by 



^(r,s) = l(V-w)(q(T,s),r), ^ 



{ q(o,s) = s. 

Along this characteristic curve, we have 



dw 

-^7(<?(t,s),t) = o, 

w(q(P,8),0) = F {a) 



which yields for r < tq 

w(q(T, s), t) = F^(s) = d q V(q(T, s), t). (2.4) 
4 



/.From (2.3)-(2.4), we obtain 



d TS q(T,s) = ±Ft>(8)d a q(T,s) - \f^{s), 
d s q(0,s) = 1. 



This yields d s q(r, s) = expf ±F* (s)t) > ifi^(a) 7^ and a g( T) s) = 1- Jritf 



if Fq(s) = when < r < tq, and then 



(s)> 



and 



V(g(T, a), r) = 2d T q(r, a) + w = 2d T q(r, M) + £ expQi^(p)r) (iftp) - itf (p)) dp + F* (s), (2.6) 

here we have used that lim(e* r (l ) H — ) = 1 — r«. 

z^o y y z> z> 

Note that q(r, M) = M such that V^g, r) satisfies the boundary condition V| g= M = 0. Together 
with (2.5)-(2.6), this yields V((?(t,«),t) = exp^(p)r) (J^(p) - i^'(p)) dp + ^(s) and g(r, a) = 

M+ [ ( exp ( -i ? o(p)r J (l - ^ ^ j + ^ ^ j dp. On the other hand, by the implicit function theorem, 

Jm\ V 2 /V F q{P) J F o\P) ) 

we can obtain the smooth function s = s(q,r) for r < tq. Therefore, V(q,r) = Fo(s(q,T)) is a smooth 

solution of (2.1) for < r < r . □ 

We now start to construct an approximate solution of (1.3) for < r = eln(l + t) < t . 
Let wo be the solution of the linear wave equation 

' dfw - Aw = 0, 
< w(0, x) = uo(x), 
dtw(0, x) = u\(x). 

Choose a C°° function x( s ) sucn th at x( s ) = 1 for s < 1 and x( s ) = for s > 2. We set, for 
< T =£ln(l + t) < T , 

u a (t, x) = e X (et)w (t, x)+ £ -{l- x (et)) X (-3eq)V(q, r). (2.7) 

By [9, Theorem 6.2.1] and Lemma 2.1, we have \Z a u a \ < C a ^£ (1 + i) _1 for r < b < r and all multi-indices 
a. Set 

J a = d?u a - (1 + u a + d t u a )Au a . 

Lemma 2.2. 

.»/« 



/ \\Z a J a \\ L 2dt <C a , b e 3 / 2 \lne\. 
Jo 



Proof. We divide the proof procedure into the following three cases. 

2 5 
Case 1. - < t < e b/e - 1. In this case, x{et) = and u a (t,x) = - x(-3£q)V(q, r). Then 
s r 



F 2 / F 2 



Ja = — - 9 9 v^ 2 v + 23 OT v/) + o 

5 



(i + *) ; 



where V(q, r) = x{—3sq)V {q, r). Since d q V has compact support, we have 

(Vd 2 V - d q Vd 2 V + 2d qT V)(q, t) 

= 9s 2 X (-3sq)x"(-Ssq)V 2 (q, r) + 27s 3 X ' \-3sq) X " \-3eq)V 2 (q, r) - 6s X '(-3eq)d T V(q, r). 

Hence \Z a J a \ < C a , b e 2 (1 + t)~ 3 + C a , 6 e 3 (1 + t)~ ip(-3sq), where is a cut-off function satisfying 

tjj(s) = 1 for s G [1, 2], and ^(s) = otherwise. 

Case 2. t < -. In this case, x( £ ^) = 1 an d « a = ei«o- This gives J a = — e 2 (w + d t w )Aw . It follows 
then from a direct computation that 

\Z a J a \ < C a e 2 (l + t)- 2 . 



1 2 

Case 5. - < t < - . A direct computation yields 



Then 



where 



u a = ew + e (1 - (r ^(-Ssg)^ - w ) 

= Jl 



Jl = ~{u a + d t u a )Au a , 

J 2 = e(8 2 - A) [(1 - X ( £ t))r- 1 x(-3 £g )(y - F )] , 
J 3 = e(d 2 - A) [(1 - x(rf))x(-3£9)(r- 1 Fo - «*)] , 
J 4 = s(d 2 - A) [(1 - x (et)) (x(-3e<?) - 1) ^o] • 

It is easy to see \Z a J x \ < C afi e 2 (1 + t)~ 2 . 

Due to (d 2 - d 2 )(V(q, r) - F (q)) = (d t - d r ){8 T V -^-) = -2& V-L- + d 2 V(-^-) 2 - d T V- 



1 + t' Tq l + t T y l + t' (l + t) 2 ' 

and V(q,0) = F (q), we have \Z a J 2 \ < C a , b e 2 \\ne\ {l+t)~ 2 . 

Moreover, by [9, Theorem 6.2.1], we have that, for any constant I > 0, if r > It, then 

\Z a (w -r- 1 F )\ < C(l + t)~ 2 . (2.8) 

On the other hand, from d t = ^ — ^r- and d r = — we obtain A = — - — (S + H)(dt — d r ) — -d r . 

t z — r z t z — r z r + 1 r 

Therefore, \Z a J 3 | < C a , b e (1 + t)~ 3 < C a , b e 2 (1 + t)~ 2 . 

Since the support of J a is in q < — — , combining the fact that, for any (f>(t, r) G C 1 , 

i + ^' 1/31 = 1 

we get the estimate |Z a J 4 | < C„, 6 e 3 (1 + i)" 1 < C Qjb e 2 (1 + t)~ 2 . 
The above analysis yields 

\Z a J a \ < C a s 2 (l + t)- 2 \\ne\. 



Collecting the estimates above, we arrive at 



\Z a J a \\ L 2 < C a , b e b ' 2 (1 + ty 1 , l<t< eWe ~ !. 
\Z a J a \\ L 2 < C a e 2 | lne| (1 + t)~ 1/2 , t < 



Consequently, 

e b ^-l 



[ \\Z a u a \\ L 2dt<C a ^ 2 \\ne\ 
Jo 



dfv - (1 + u + d t u)Av = -J a + (v + d t v)/\u a , 



(2.11) 



and Lemma 2.2 is proved. □ 

For later reference, we cite a result that was shown in [17]. 

Lemma 2.3. Let f(t,x) S C 1 (M 4 ) only depend on (t,r). Moreover, assume that supp/ C {(t,x) : r < 
M + t}. Then 

||(l + |*- r |)-7|| L2 <C1|d r /|U 2 . 

Based on these preparations, we next establish: 

Proposition 2.4. For sufficiently small e > and < r = eln(l + 1) < b < To, (1.3) has a C°° solution 
u(t,x) which satisfies, for all \a\ < 2, 

\Z a d(u - u a )\ < C b e 3 ' 2 i lne|(l + t)-\l + \t - H)" 1 / 2 . (2.10) 
Proof. Set v = u — u a - Then one has 

r d 2 v - 

1 v(0,x) = d t v{0,x) = 0. 
We make the induction hypothesis that, for some T < e h l e — 1, 

\Z a dv\<e{l+t)- 1 {l + \t-r\)- 1/2 , \a\ < 2, i < T, (2.12) 
which then further implies that, for |a| < 2 and t < T, 

\Z a v\ < Ce(l + + |t - r|) 1/2 . (2.13) 

To verify the validity of (2.12), we will prove that, for sufficiently small e > 0, 

\Z a dv\< S -{l+ty\l + \t-r\)- 1 l 2 , |a|<2, t < T, (2.14) 

and apply the continuity method to obtain eln(l + T) = b. 

Applying Z a to both sides of (2.11) and using [Z a ,d 2 - A] = ^ C^Z (d 2 - A) yields, for any 

I/3KM 

H < 4, 

(a t 2 -A)Z« v = Z«G- ^ C Q/3 ^G, (2.15) 

l/3|<l«l 

where 

G = (u + d t u) Av - J a + (v + d t v) Au a . 
7 



Since 

Z a G = {u + d t u)Z a Av+ Z ai (u + d t u)Z a2 Av + Z a [-J a + (v + d t v)Au a } 



|ai|>l 



we have from (2.15) that 

(<9 4 2 - (l + u + d t u)A) Z a v = F, (2.16) 

where 

F= Z ai {u + d t u)Z a2 Av + Z a [-J a + (v + d t v)Au a }- ^ C af} Z p G + (u + d t u) [Z a , A] v. 

a 1 +a 2 =a \P\<\<*\ 
|ai|>l 

Next we derive an estimate of ||9Z a v||i2 from equation (2.16). Define the energy 
£(t) = JV ( (\d t Z a v\ 2 + (l + u + dtu)\VZ a v\ 2 )dx. 

2 M<4 7R3 

Multiplying both sides of (2.16) by dfZ a v {\a\ < 4) and integrating by parts, and noting that \d@u\ = 
\d^u a + d^v\ < Cb£ (1 + t) _1 (|/3| = 1, 2), which follows from the construction of u a and assumption (2.12), 
we arrive at 

E'(t)<^-E(t)+J2 I \F\\d t Z°v\dx. (2.17) 
We now treat each term arising in the integration of V / \F\ | W dx separately. 

\a\<4 J *' 

{A) The term V / \Z ai (u + d t u)Z a2 Av\ \d t Z a v\ dx. First, we note that: 

l«ll>l 

(i) By assumption (2.13), for |/?| < 2, we have 

| ( l + | t _ r |)-i Z P (Uo + u) |<^L. (2.18) 

(ii) By (2.18) and (2.9), for \a x \ + \a 2 \ = \a\ < 4 with > 1, we have 
/ \Z ai {u + d t u)(Z a2 Av)d t Z a v\dx 

JR 3 

<C b f \(Z ai u a + Z ai d t u a )(Z a2 Av)dtZ a v\dx+ f \(Z ai v + Z ai d t v)(Z a2 Av)d t Z a v\dx 
Jr 3 Jr 3 

<3L E (t) + C V / \{l + \t-r\)- 1 Z a ^v\\Z'<dv\\d t Z a v\dx 

|7|<|a 2 | + l 

+ C b ^ / |^ ai 9t«| \Z-*dv\ \d t Z a v\dx. (2.20) 

|7|<I«2| + 1^ R3 

Note further that there is at most one number larger than 2 between \ai\ and |7|. If \ai\ > 2, then |— y | < 2 
and, by Lemma 2.3 applied to (1 + \t — r\)~ 1 Z ai v and by assumption (2.12), we arrive at 

/ |(1 + \t-r\)- 1 Z a 'v\ \Z^dv\ \d t Z a v\dx+ f \Z a ^d t v\\Z^dv\\d t Z a v\dx<^-E(t). (2.21) 

JR 3 ' JR 3 1+t 



If |7| > 2, then |ai| < 2 and it follows from (2.13) that |(1 + \t - r\)~ 1 Z a ^v\ < Ce(l + £) - 1 (1 + 1 1 — r | ) - 1 / 2 
which leads to 

/ \{l + \t-r\)- 1 Z ai v\\Z'<dv\\d t Z a v\dx+ [ \Z ai d t v\ \Z^dv\ \d t Z a v\dx < E(t). (2.22) 
jr s Jr s 1 +t 

Inserting (2.21)-(2.22) into (2.20) yields 

V f \Z ai (u + d t u)Z a2 Av\\d t Z a v\dx < ^-E(t). (2.23) 

JR3 l + t 

ai+a2-cv 
l«i|>l 

(B) The terms J R3 \Z^{(u + dtu)Av) ■ dtZ a v\dx with \/3\ < \a\. We only need to treat the term 
Jks \ (u + d t u)Z^Av ■ d t Z a v\dx, since the other terms have been estimated in (A). 

By (2.12), we have 

f \(u + d t u)Z Av d t Z a v\dx 

^ C zZ f \{± + \t-r\r 1 u\\Z~<dv\\dtZ a v\dx+ [ \d t u\\Z l3 Av\\d t Z a v\dx 

<^f- t E(t). (2.24) 

(C) The terms f R3 \d t Z a v\ \Z /3 J a \ dx with |/3| < \a\ < 4. In this case, we have 

/ \dtZ a v\\zPj a \dx< \\ZPj a \\ L ,y/E(tj. (2.25) 

JR* 

(D) The terms / R3 \Z?((v + dtv)Au a )\ \dtZ a v\ dx with |/3| < \a\ < 4. A direct computation yields 
/ \zP((v + d t v)Au a )\\d t Z a v\dx 

JRS 

<C V / (\{l + \t-r\)-\Z^v){Z^du a )\ + \{Z^d t v){Z^du a )\)\d t Z«i 

. JVL 3 



*v\ dx 



\Pi 1 + 1/32 |<|/3| + 1 
l/3i|<l/3| 



<^f- t E{t). (2.26) 
(E) The term J R3 \(u + d t u) [Z a , A] v\ \d t Z a v\ dx. Since [Z a , A] = ^ C al3 d 2 Z^ we have 

\P\<\a\-l 

f \(u + d t u)[Z a ,A}v\\d t Z a v\dx <C V / \u + d t u\\d 2 Z l3 v\\dtZ a v\dx 
Jrs \p\<\a\-i J *' 

^ C Y, I \(l + \t-r\)- 1 (u + d t u)\\dZ''v\\d t Z a v\dx 
Jrs 

\0\<\a\-l 
|7|<I/3| + 1 

<^- t E(t). (2.27) 



Substituting (2.23)-(2.27) into (2.17) yields 

C b e 



E'(t) < T ±-E(t)+ ]T \\ZPj a \\ L ,y/E(t). (2.28) 

l/3|<4 

Thus, by Lemma 2.2 and Gronwall's inequality we obtain 

\\dZ a v\\ L 2 < C 6 e 3 / 2 |ln £ |, |a|<4, 

and then 

\\Z a dv\\ L 2 < C b e 3/2 |lne|, \a\ < 4. (2.29) 
By (2.29) and the Klainerman-Sobolev inequality (see [9] or [14]), we have 

\Z a dv\ <C b £ Z/2 \\n£\(l + t)- 1 {l + \t-r\)- 1 /' 2 , |a|<2, t<T. (2.30) 

This implies that, for e > small enough, 

\Z a dv\<^(l + t)- 1 (l + \t-r\)- 1 / 2 , \a\<2, t<T. 

Therefore, we have completed the proof of (2.12) and, together with (2.30), the proof of (2.10). □ 
Proposition 2.4 immediately gives lim s ln(l + T s ) > tq and, therefore, 

UmelnT e >T . (2.31) 

Remark 2.1. The analysis of this section can be directly applied to problem (1.6) with general initial data 
so that a lower bound on the lifespan T £ is obtained as in [5] . 

§3. The upper bound on the lifespan T e and proof of Theorem 1.1. 

In this section, we focus on the upper bound on T s . Here some ideas are inspired by [10] and [6]. Since 
equation (1.3) contains the solution u and its derivatives simultaneously rather than the derivatives of u 
only, as in [10], and the function u only, as in [6], respectively, we have to be more careful in computations 
and also need to estimate more quantities. Thanks to the estimate of Z a (u — u a ) with |a| < 2 in (2.13), 
we observe that \Z a (u — u a )\ < Cs (1 + 1)~ holds near the light cone which plays an important role in the 
analysis below. 

Let U = ru and c 2 (u) = 1 + u + dtu. Due to radial symmetry of u, (1.3) can be rewritten as 



r d 2 u - c 2 d 2 u = o, 

\ U(0,r) = eru , d t U(0, r) = eru\. 

Define two operators 

L\ = d t + cd r , L 2 = d t - cd r . 

We also set 

wi = L2d r U = (d t — cd r ) d r U, w 2 = L\d r U = (d t + cd r ) d r U, 

which yields Ot r U = and a r v = — • 

Note that 

LiL 2 = d 2 - c 2 d 2 - (L lC )d r , L 2 L X = d 2 - c 2 d 2 . + [L 2 c)d T . 
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(3.1) 



Then 

Liw\ = LiL 2 d r U = —- — w\ + - — (w 2 + d t u + -L 2 u) - - — (d t u + -L 2 u), (3.2) 
Arc Arc c Arc c 

1 ItJ^ f UU\ f 

L 2 w 2 = L 2 Lid r U = - — Wo + - — {—w\ - d t u + -Liu) + - — (d t u L\u). (3.3) 

Arc Arc c ire c 

Due to d r c = c'd r u = —d r u — - — d t u + - — (w\ + w 2 ), we have 
2c 2cr Arc 

L\W\ + wid r c = —^-(2w 2 — dtu H — L\u) — ^-(dtu H — L 2 u), 
Arc c Arc c 

L 2 w 2 - w 2 d r c = ^-(—2wi + d t u + -L 2 u) + ^-<d t u Liu) 

Arc c Arc c 

and 

d(\wi\(dr — cdt)) = sgn w\ {L\W\ + w\d r c) dt A dr 

r W\ r w 2 r i 

= sgn w\ (2w 2 — d t u + -Liu) — - — (d t u + -L 2 u) dt A dr, (3.4) 

I Arc c Arc c J 

d(\w 2 \(dr + cdt)) = sgnw 2 (L 2 w 2 — w 2 d r c) dt A dr 

r w 2 r wi r i 
= sgn w 2 - — (— 2wi + dtu H — L 2 u) + - — (dfU L\u) dt A dr. (3-5) 



t 

A- 




Figure 1. 



/.From §2 we have that when eln(l + X&) = b < to and b > is a fixed constant, (1.3) has a C°° solution 
for t < Tb. Choosing e sufficiently small such that 1/e < e b l £ — 1, we define the characteristics by 
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dr/dt = ±c, respectively, passing through (A, 0) in the (r, t) plane. Let D be the domain which is bounded 

2 F"(o ) 

by T+ and T+ , (see Figure 1), where p is chosen such that r(p ) = — — - In — - — 2 — <L- — = Tq or 
2 

— - — - = t . (Another possibility is that there exists no point po G (— M, M) such that r(po) = i~ . In this 
F o (Po) 

case, by a minor modification which uses the fact that, for any fixed Sq > sufficiently small, there exists 
Po € (— M, M) such that t(pq) < to + So, an analogous proof still works). Obviously, is the straight 
line r = t + M. 

Fix a positive constant R satisfying R > To, and assume that u is a solution of (1.3) (at least) for 
t<T < exp(R/e). We define 

A(t) = sup / \wi (s, r)\ dr, 

l/e<s<t J(s,r)eD 

B(t) = sup s\d s u(s,r)\, 

l/e<s<t 
{s,r)€D 

C(t)= sup s 3 / 2 \w 2 ( Sl r)\, 

l/e<s<t 

D(t)= sup s 3/2 \L 1 u(s,r)\. 

l/e<s<t 
(s,r)eD 

Lemma 3.1. There exists a constant E > such that, for e > sufficiently small, 

771 

A(l/e)<—, B(l/e)<Ee, C(l/e)<E 2 e 2 , D(l/e) < Ee. (3.6) 



Proof. The equation r = r(t) of for AG [po — 1, M] is 

f dr(t) 



dt 

r(0) = A. 



(u(t,r(t))) = c(t), 



This gives 

r(t) - A = / (c(s) - l)ds + t. 



Because of |c(i)-l| < C 6 (|u(t, r(t))| + \d t u(t,r(t))\) < C fe £r(l+t)- 1 (l + |*— r(t)l) 1 / 2 for < r = £ln(l+i) < 
b < t , 

|r(t) - t\ < I A| + C b [ s(l + sy 1 (1 + \s - r(s)l) 1 / 2 
Jo 

<m + b + C b £(l + s) _1 |s-r(s)|(is, (3.7) 



where mo = max{|po — 1 1 , M}. It follows from Gronwall's inequality that 

\r(t)-t\<C b . 

Therefore, 



\t + M-r(t)\<C b , (3.8) 
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which means the distance between r^" o _ 1 and is bounded. 
The equation r = r(t) of T~ is 

\ r(0) = fi, 

and then \r(t) + 1 — fi\ < C b can analogously be shown. 

For (t,r), (£',/) G T" n £> for ^ G R and (f , r) G 17^, (i',r') G rj, we obtain 

\t - t'\ < ^ [\t + r - H + \t' + r' - m| + \t - r - A| + \t' - r' - X'\ + |A - A'|] < C b . (3.9) 

We now start to prove (3.6). Since Wi = d t u + rd tr u — cd r u — crd^u, one has, for t < e b / £ — 1, 

|wi(t, r)\ < C b e (l + \t- r|)" 3/2 . (3.10) 

By (3.8)-(3.10), we thus arrive at 

/ \wi(s,r)\dr < C b e. (3.11) 

In addition, because of \d t u{t, r)\ < C b s (1 + t)' 1 (1 + \t - r])" 1 / 2 for t < e 6 / 6 " 1 , we have, for (1/e, r) G D, 

< C b s. (3.12) 



c7, ti[ -. i 



Note that 

W2 (t, r) = (d t + cd r )u + —(S + H)d r u + (c - l)rdlu 
r + t 



which implies \u)2{t,r)\ < C b e(l + t) 1 for (t.r) G D and sln(l + t) < b in view of (2.13). Together with 
(3.3), this yields 

\L 2 w 2 \<0-^. (3.13) 
Because of w 2 (t, t + M) = and (3.9), we obtain from (3.13) that 

C(l/e) < C b e 5 / 2 < C 6 e 2 . (3.14) 

Finally, from L x u = — !— (S + H)u + (c - l)d r u we have that |Li«(t, r)\ < C b e(l + t)~ 2 for t < e h ' e - 1 

T -\- t 

and, therefore, 

D{l/e) < C b e 3 ' 2 < C b e. 
Collecting (3.11)-(3.12) and (3.14) completes the proof of (3.6). □ 

Based on Lemma 3.1, we will use the continuity method to establish the upper bound on T e . To this 
end, we assume that, for < t < T' < T, 



A(t) < Ee, B(t) < 2Ee, C(t) < 3E 2 e 2 , D{t) < 2Ee. 
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(3.15) 



We now start to verify: 

Lemma 3.2. Under the assumption (3.15) and if e > is sufficiently small, then we have, forl/e < t < T' , 



A{t) < \Ee, B(t) < Ee, C(t) < ^E 2 e 2 , D(t) < Ee. 



(3.16) 



Proof. We first estimate u(t, r) for 1/e < t < T' and (£, r) £ D. Let (£, r) £ T~^. We integrate L\u along 
from time 1/e to t. From assumption (3.15), it is readily seen that |u(£,r)| < Ce 3 ! 2 . Hence the definition 
c(t, r) = Vl + u + dtu is possible and assumes a value close to 1. Therefore, t \d r u(t, r)\ < Ce which yields 
\u(t,r)\ < Ce(l (1 + \t- r\) for 1/e < t < T and (t,r) e D. 

Similar to the proof (3.8) and (3.9) we also have \r — t\ < C and \t — t'\ < C, where t and t' are same as 
in (3.9). 

Now we estimate A(t). For 1/e < t < T", by equation (3.4) and Green's formula, we have 



{t,r)£D 



\wi(t,r)\dr 



< / \wi(l/e, r)\dr+ // ( 2w 2 - d t u + -L x u ) - ( d t u + -L 2 u ) (s.r)dsdr 

J(i/e,r)eD J J 4rc V c ) Arc \ c ) 



< 



'(l/e,r)6 

Ee 



II Lrc( 



l/e<s<t 
(s,r)£D 

r 

2w 2 — dtu H — L\u 



)- — ( 

/ Arc \ 



d t u + -L 2 u 

c 



) (s,r) 



dsdr. 



(3.17) 



l/e<s<t 
(s,r)GD 



2Ee 3E 2 e 2 2Ee 

By the induction hypothesis (3.15), we have \d s u(s,r)\ < , Iw2(s,r)| < — tttt, and \Liu(s,r)\ < -^-nr 

S S ' s ' 

5Ee 

for 1/e < s < T' and (s,r) G D. This immediately gives \L 2 u{s,r)\ < . Note also that \r — s\ < C 

holds for s > 1/e. We then have r > s/2 and 

2 [dtu H — L 2 uj (s, r) dsdr < J J 



IS 



Arc 



E 2 e 2 (2Ee 6Eer\ 
s 3/2 r y s S J 



dsdr 



l/e<s<t 



< 25£ 3 e 3 



f — f 

Jl/e « 5/2 J(s,r)eD ' 



l/e<s<t 
(s,r)eD 

dr < CE 3 e 9/2 < 



Ee 
12' 



(3.18) 



where the generic constant C > is independent of e. 
Similarly, 



w\ r 

(2w 2 — dtu H — Liu) I (s, r)dsdr < 

Arc c 



K(s,r)| (6E 2 e 



3r 



,3/2 



2Ee . 2Eer 
s 



S 3/2 C 



dsdr 



l/e<s<t 
(s,r)eD 



l/e<s<i 
(s,r)Efl 



<AEe /' 4t7 / |wi(a,r)|rdr <CE 2 e 5/2 . 

Jl/e 8 1 J{s,r)eD 



(3.19) 



Substituting (3.18)-(3.19) into (3.17) yields 



A(t)<^ + ^ + CE 2 eV 2 , 
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which implies that A(i) < | Ee for e > sufficiently small. 
Next we estimate B(t). Note that u satisfies the equation 



L 2 d t u = -- Wl + - d t u. (3.21) 
r r 




Figure 2 



We integrate equation (3.21) along the characteristics r which insects at the point (t\r'). For 
t' > 1/e, we denote by D\ the domain bounded by T~, the line {t = £}, and T^- (see the Figure 2) which 
are denoted by r 1; T 2 and r 3 , respectively. We have 

J J sgnu/i - d t u + ^Liuj - (d t u + ^L 2 uj^ (s, r) dsdr = (^J + J ^ \wi\ (dr - cdt) 

which implies 

/ |wi| (dr — cds) 

J J (21V2 — dtu + -Liuj — (dtu + -^2uj (s,r) dsdr + J \wi \ (dr — cdt) 

W\ / r \ W2 / r \ I f 

- — 2w 2 - d t u + -Liu - - — d t u + -L 2 u (s,r) dsdr + / \w 1 (t,r)\dr 
Arc V c ) ire \ c )\ i(t, r ) eDl 



(s,r)SDi 



< 



l/e<s<t 
(s,r)€D 

1 ^ 2 „ 5 „ 
< - + - £e = - Se. 
~ 6 3 6 
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This yields 



\wi(s,r(s))\ 



, f KM«))| 2 /r 



« 7r 



(dr — cds) 



< — Ee. 
~ 6t 



s 1 

Note that r(s) > - > — holds. ^,From (3.21) we then have 

z AS 



\d t u{t,r)\< [ 
Jt> 



\cwi\(s,r(s)) /"* \cd t u\(s,r(s)) 



r{s) 



-ds + 



r(s) 



■ ds 
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<—Ee + 3e I \d t u(s,r(s)\ ds. 



(3.22) 



By Gronwall's inequality, we obtain 



1 1 Ee 
\u(t,r(t))\<— Eee Cs <^. 



For t' < 1/e, we have t <t' + \t-t'\ < 2/e, and Section 2 tells us \td t u(t, r)\ < Ee. Therefore, B(t) < Ee 
holds. 

Finally, we estimate C(t). We write (3.3) as 



L 2 w 2 = aw 2 + b, 



(3.23) 



where 



Arc 



w 2 - Wi 



dtu H — Liu) , b = wi (dtu — -Liu) . 

c / Arc V c / 



Integrating (3.3) along r as above, we have 



\w 2 (t,r)\< / \au)2 + b\(s,r(s))ds. 



(3.24) 



2i 

Noting f > t — \t — t'\ > — , we obtain 

3 



j\b(s,r(s))\ds< £ 



1 



- — wid t u 
Arc 



(s, r(s)) ds + 
2Ee f \ Wl {s,r{s))\ ^ ^ 5E 2 e 2 



(s, r(s)) ds 



2Ee f , 

*V2 J t , r(s) aS S 3t 3 /2 



(3.25) 



and 



i 



aw 2 \ (s, r(s))\ ds < 



E 2 e 2 /•* \w 2 -w 1 -d t u+ ^L lU \(s,r(s)) CE 3 e i 
' ds < 



V' 2 j t 

Substituting (3.25)-(3.26) into (3.24) yields 



r{s) 



t 3/2 



\w 2 (t 7 r)\<J^ \aw 2 + b\ (s,r(s)) ds < 



2J2 



5E 2 e 
2*3/2 
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(3.26) 



5E 2 e 2 

and C(t) < — - — is shown. 
Eventually, we have that 

L2L1U = d r u + —(L2u)(d r u) + — ( — wi H — dtu) d r u. 

r 2c 2c V r r ) 

We integrate this equality along T~ from t' to t and obtain 

. T . [ f 6Ee , f* 15E 2 e 2 , f* f \ Wl (s,r(s))\ Ee \ 3Ee , 

CEe Ee 
~ ~ W 2 

and, therefore, D(t) < Ee. □ 

Based on Lemmas 3.1-3.2, we will use the following blowup result to establish the upper bound on the 
lifespan T e . 

Lemma 3.3. Let w be a solution in (0,T) of the ordinary differential equation 

dw . . 2 

— = a {t)w + a\(t)w + a 2 (i) 
with aj continuous on [0, T] and a > 0. Let 



*3 

rT \ / t-T 



\a\(t)dt . 



K= \L |a2(t)l ^J exp yj 

lfw(0) > K, then 

(j^ ao(t) dtj exp ^- \ ai (t)\dtj < (w(O)-iTT 1 . 

Proof. The proof can be found in [8, Lemma 1.4.1] or [9, Lemma 1.3.2]. □ 
Next we show that 

TmTslnTe < t q . (3.27) 

e— >0 

Proof of (3.27). Along the characteristics r+ , wi(t,r(t)) satisfies 

dwi , r 1 o u>i ( „ r r \ w 2 /„ 

— (t, r(t)) = = -— Wl + — (™ 2 + ft« + -L aU ) - — (d t u + -L 2 u) . 

We fix < fi < 1 and set t e = e^ T "/ e - 1 and w(t,r(t)) = r(i)) exp^- ^(s,r(s)) ds^j . Then 



-(t, r(t)) = a (i)w 2 + oi(t)w + o 2 (t), 



where 



aoW = 47(^ exp (I^ (s ' r(s))ds )' 



w 2 + d t u 
fl iW = 4rc 



02 (*) = + ^ L 2w) exp^-^ ^(s,j"(s))ds^(t,r(t)). 
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By (3.15), we have that, for < t < T, 



. Ee 
m\ < -nr, ml < 



t 2 > f 5/2 

which yields 

/ |ai|ds<£e 2 , / |a 2 |ds < CE 2 s 9/2 . 

Jtr Jt r 



This further yields 

if 



\a 2 {t)\ d?j exp ^ \ ai {t)dtj =0(e 9/2 ). (3.28) 



5 

Duetoi e > 2/e and |3e(£ e — r(£ e )| < 1, by the definition of u a in (2.7), we have that ii a (t e ) = — ; — rV(t E ) 
V(r(t £ ) — t £ , jUTo) holds. In addition, it follows from Proposition 2.4 that 



r(t e ) 

|Z> - «„)| < C a ^' 2 | lne| (1 + ty 1 (l + \t- r\) 1/2 . 

Therefore, 

wi(t £ ) = (d t u + rd tr u)(t £ ) - [c(2d r u + rd 2 u)] (t £ ) 

= (d t u a + rd tr u a ){t e ) - [c(2d r u a + rd 2 u a )] (t e ) + o(e) 

= -2ed 2 V(t £ ) + ^-d qT V{t £ ) + o(e) 

= -2ed 2 q V(t £ ) + o(e) 

which yields w(t £ ) = 2ed 2 V{t £ ) + o(e). 
By Lemma 3.3, we obtain 



a (t) dtj exp^- jf |ai(t)|dt) < (w(t e ) - -K") \ 



and then 

(LwC & <s ' *»*) *) eip (- r |ai<t)| *) * <ffi(t - ) - ^ ■ 

It follows 

( In(T e + 1) - ln(e^'>/ e + e/2)) (1 + 0(e)) < 4{2ed 2 V{t £ ) + 0(e))" 1 exp(£ £ 2 ). 
In order to obtain (3.27), we need to estimate V(t £ ). Set q(r) = r(t) — t. Then 
( dq t + l 



dr r(t)(l + c(t)) 
9(0) = Po- 



(V-d q V)(q,T) + 0(e^ 2 \hiB\), 



Recalling that, in Lemma 2.1, q(r, po) is the characteristics of the approximate equation which satisfies 
q(0, po) = po, and then comparing q with q yields 

d{q-q)_( 2(t+l) ,U-9/,_ . , (V-d q V)(q,r) (V - d q V)(q, r) 



dr 



+ o(e 1/2 |lne| 



2 2 



I g(0)-g(0,p ) = 0. 
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By Gronwall's inequality, we have q(r) — q{r,po) = O (e 1 ' 2 ! lneQ when r < R. Thus, d 2 V(t e ) = 
dqV(q(nTo, po), po) + (s 1 ^ 2 ] In e | ) . From the expression for V in Lemma 2.1, we have that dgV(q(r, s), r) = 
F"(s) F"(o ) / \ 

- °, \ , and consequently d?V(t e ) = r ° — . + O ( e 1/2 \ lnel ) . Moreover, one has F{!(po) > by 
O s q{T,s) q d s q{fiT ,p ) \ J 

Remark 1.2. Therefore, 

2d s q(p,T , pa) 



and as — > 1, 



limeln(T e + l)-^T < 
e^o Fo(Po) 



hmelnTe < t H =77 — r — = r . 



This finishes the proof of (3.27). □ 

Proof of Theorem 1.1. Follows now directly from (2.31) and (3.27). □ 



References 

1. S. Alinhac, Blow up of small data solutions for a class of quasilinear wave equations in two space dimensions, Ann. of 
Math. (2) 149 (1999), 97-127. 

2. , Blow up of small data solutions for a class of quasilinear wave equations in two space dimensions. II, Acta Math. 

182 (1999), 1-23. 

3. , An example of blowup at infinity for quasilinear wave equations, Asterisque 284 (2003), 1—91. 

4. D. Christodoulou, Global solutions of nonlinear hyperbolic equations for small initial data, Comm. Pure Appl. Math. 39 
(1986), 267-282. 

5. Ding Bingbing, I. Witt, and Yin Huicheng, On the lifespan of and the blowup mechanism for smooth solutions to a class 
of 2-D nonlinear wave equations with small initial data (2012), arXiv:1210.7980. 

6. Ding Bingbing and Yin Huicheng, On the blowup of classical solutions to the 3-D pressure-gradient systems, J. Differential 
Equations 252 (2012), 3608-3629. 

7. P. Godin, Lifespan of solutions of semilinear wave equations in two space dimensions, Comm. Partial Differential Equa- 
tions 18 (1993), 895-916. 

8. L. Hormander, The lifespan of classical solutions of nonlinear hyperbolic equations; Mittag-Lefner report no. 5 (1985). 

9. , Lectures on nonlinear hyperbolic equations, Math. Appl., vol. 26, Springer, Berlin, 1997. 

10. F. John, Blow-up of radial solutions of utt = c 2 (ut)Au in three space dimensions, Mat. Apl. Comput. 4 (1985), 3-18. 

11. F. John and S. Klainerman, Almost global existence to nonlinear wave equations in three space dimensions, Comm. Pure 
Appl. Math. 37 (1984), 443-455. 

12. M. Keel, H. Smith, and C. D. Sogge, Almost global existence for quasilinear wave equations in three space dimensions, 
J. Amer. Math. Soc. 17 (2004), 109-153. 

13. S. Klainerman, The null condition and global existence to nonlinear wave equations, Nonlinear systems of partial differ- 
ential equations in applied mathematics, Part 1 (Santa Fe, N.M., 1984), pp. 293-326; Lectures in Appl. Math., vol. 23, 
Amer. Math. Soc, Providence, RI, 1986. 

14. , Remarks on the global Sobolev inequalities in the Minkowski space M n ^ 1 , Comm. Pure Appl. Math. 40 (1987), 

111-117. 

15. S. Klainerman and G. Ponce, Global, small amplitude solutions to nonlinear evolution equations, Comm. Pure Appl. 
Math. 36 (1983), 133-141. 

16. Li Ta-tsien and Chen Yun-mei, Initial value problems for nonlinear wave equations, Comm. Partial Differential Equations 
13 (1988), 383-422. 

17. H. Lindblad, On the lifespan of solutions of nonlinear wave equations with small initial data, Comm. Pure Appl. Math. 
43, no. 4, 445-472 (1990). 

18. H. Lindblad, Global solutions of nonlinear wave equations, Comm. Pure Appl. Math., 45 (1992), 1063-1096. 

19. , Global solutions of quasilinear wave equations, Amer. J. Math. 130 (2008), 115-157. 

20. H. Lindblad, M. Nakamura, and C. D. Sogge, Remarks on global solutions for nonlinear wave equations under the standard 
null conditions, J. Differential Equations 254 (2013), 1396-1436. 

21. J. Metcalfe and C. D. Sogge, Global existence for Dirichlet-wave equations with quadratic nonlinearities in high dimen- 
sions, Math. Ann. 336 (2006), 391-420. 

22. , Global existence of null-form wave equations in exterior domains, Math. Z. 256 (2007), 521—549. 

23. Wu Sijue, Almost global wellposedness of the 2-D full water wave problem, Invent. Math. 177 (2009), 45-135. 

24. , Global wellposedness of the 3-D full water wave problem, Invent. Math. 184 (2011), 125-220. 

25. Yin Huicheng, The blowup mechanism of small data solution for the quasilinear wave equations in three space dimensions, 
Acta. Math. Sinica, English Series 18 (2001), 35-76. 

19 



26. , Formation and construction of a shock wave for 3-D compressible Euler equations with the spherical initial data, 

Nagoya Math. J. 175 (2004), 125-164. 

27. Zhou Yi and Xu Wei, Almost global existence for quasilinear wave equations with inhomogeneous terms in 3D, Forum 
Math. 23 (2011), 1113-1134. 

Department of Mathematics and IMS, Nanjing University, Nanjing 210093, P.R. China 
E-mail address: balancenjustayahoo.com.cn 

Mathematical Institute, University of Gottingen, Bunsenstr. 3-5, D-37073 Gottingen, Germany 
E-mail address: iwittOuni-math.gwdg.de 

Department of Mathematics and IMS, Nanjing University, Nanjing 210093, P.R. China 
E-mail address: huicheng@nju.edu.cn 



20 



